FLOW OF ELECTRORHEOLOGICAL FLUIDS UNDER THE 
CONDITIONS OF SLIP ON THE BOUNDARY. 
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Abstract. We derive general conditions of slip of a fluid on the boundary. Under these 
conditions the velocity of the fluid on the immovable boundary is a function of the normal 
and tangential components of the force acting on the surface of the fluid. A problem on 
stationary flow of an electrorheological fluid in which the terms of slip are specified on one 
part of the boundary and surface forces are given on the other is formulated and studied. 
Existence of a solution of this problem is proved by using the methods of penalty functions, 
' monotonicity and compactness. It is shown that the method of penalty functions and the 

, Galerkin approximations can be used for the approximate solution of the problem under 

■ consideration. 

m 



1. Introduction 

O . 

, Electrorheological fluids are smart materials which are concentrated suspensions of po- 

CN ' larizable particles in a nonconducting dielectric liquid. In moderately large electric fields, 

■ the particles form chains along the field lines, and these chains then aggregate to form 

I columns [9]. These chainlike and columnar structures cause dramatic changes in the rheo- 



logical properties of the suspensions. The fluids become anisotropic, the apparent viscosity 



Q^. (the resistance to flow) in the direction orthogonal to the direction of electric field abruptly 



increases, while the apparent viscosity in the direction of the electric field changes not so 
drastically. 

The chainlike and columnar structures are destroyed under the action of large stresses, 
and then the apparent viscosity of the fluid decreases and the fluid becomes less anisotropic. 

■ On the basis of experimental results, the following constitutive equation was developed 
: in [3]: 

■ aij{p,u,E) = -p5ij + 2ip{I{u),\E\,fi{u,E))eij{u), i,j = l,...,n, n = 2 or 3. 

(1.1) 

Here, aij{p,u, E) are the components of the stress tensor which depend on the pressure p, 
the velocity vector u = {ui, . . . , Un) and the electric field strength E = (Ei, . . . , En), 6ij are 
the components of the unit tensor (the Kronecker delta), and £ij{u) are the components of 
the rate of strain tensor 

Moreover, I{u) is the second invariant of the rate of strain tensor 

n 

I{u)=Y.{e.,{u))\ (1.3) 
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and is the viscosity function depending on I{u), \E\ and fi{u, E), where 

(a(u E))(x) = = irl=,u^i-)E^{x)r 

So fi{u, E) is the square of the scalar product of the unit vectors and The function /x 
is defined by (1.4) in the case of an immovable frame of reference. If the frame of reference 
moves uniformly with a constant velocity u = (?2i, . . . , ti„), then we set: 

As the scalar product of two vectors is independent of the frame of reference, the constitutive 
equation (1.1) is invariant with respect to the group of Galilei transformations of the frame 
of reference that are represented as a product of time-independent translations, rotations 
and uniform motions. 

The presence of the function n in the constitutive equation (1.1) is connected with the 
anisotropy of the electrorheological fluid under which the viscosity of the fluid depends on 
the angle between the vectors of velocity and the vector of electric field. 

The function \x defined by (1.4), (1.5) is not specified at £^ = and at it = 0, and there 
does not exist an extension by continuity to the values of u = and £' = 0. However, 
at -B = there is no influence of the electric field, and the function [i{u^ E) need not to 
be specified at £^ = 0. Likewise, in case that the measure of the set of points x at which 
u{x) = is zero, the function /x need not also be specified at u = 0. But in the general the 
function fi can be defined as follows: 

ai + u{x) + u E{x) \2 
^a^/n + \u{x) +u\'' \E{x) 

where / denotes a vector with components equal to one, and a is a small positive constant. 
If u{x) 7^ almost everywhere in Q,, one can choose a = 0. 

The viscosity function is identified by approximation of fiow curves, see [3], and it was 
shown in [3] that it can be represented as follows: 

^{I{u), \El^l{u, E)) = h{\Elii{u, E)){X + /(tx))-5 + i;{I{u), \E\, ii{u, E)), 

(1.7) 

where A is a small parameter, A > 0. 

The equations for the functions E and {p,v) are separated, (see [3]). Because of this, we 
assume here and thereafter that the function of electric field E is known. 

Various problems on stationary flow of electrorheological fluids under mixed boundary 
conditions such that velocities and surface forces are prescribed on different parts of the 
boundary are investigated in [3] . This formulation assumes that the fluid adheres to a hard 
boundary, that is the velocity of the fluid on the hard boundary is equal to the velocity of 
the boundary. 

But at some conditions, wall effects appear, the velocity of a fluid on the hard boundary 
can be different from the velocity of the hard boundary. In particular, hard particles of 
electrorheological suspensions may slip along the hard boundary. 

It was shown experimentally that magnetic suspensions, whose conduct is similar to the 
conduct of electrorheological fluids, exhibit wall effect, see [5]. This effect depends both on 
the surface roughness of the wall and on the force pressing the particles against the surface 
of the wall. 



^, , , / ai + u(x) + u E(x) \2 
Li(u,E)(x) = i — ^ — -Vv^ 1^ M , 1-6) 
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In Section 2, we derive the boundary conditions of slip. In Sections 3 and 4, we formulate 
a boundary value problem on stationary flow of the electrorheological fluid under the con- 
dition of slip on the boundary and present a theorem on the existence of a solution of this 
problem. Sections 5 and 6 are devoted to the proof of the existence result and construction 
of approximate solutions by using the method of penalty functions. In Section 7, we show 
that Galcrkin approximations can be used for approximate solution of our problem. 

Since the constitutive equations of nonlinear viscous and viscous fluids are partial cases 
of the equation (1.1), the results presented in this paper can be applied to nonlinear viscous 
and viscous fluids. 

2. Frictional force and the velocity of slip on a hard boundary. 

Let Q C be a domain in which a fluid flows. Let S be the boundary of 0, and Si be 
a part of S which corresponds to a hard immovable wall. We assume that the fluid slips on 
^i. Let F{s) = -Pi(s)Ci be an external surface force acting on the fluid. Here Q are 

unit vectors directed along the coordinate axes Xi, Fi scalar functions of points s of iSi. 

We represent the function F in the form 

F{s) = F''{s)+F-{s), seSi, (2.1) 
where and F'^ are the normal and the tangential surface forces. 

n 

F'^is) = F.isMs), F,{s) = J2 F^{sMs), (2.2) 

i=l 

n 

F^{s) = Fis) - F-'is) = FriisXi, Friis) = Fi{s) - F,{sMs), (2.3) 

i=l 

where i' = (ui, . . . ,Un) is the unit outward normal to Si. 

Analogously, the velocity vector u on the boundary is represented in the form 

n 
i=l 

n 
i=l 

n 

u^s) = u{s) - u'is) = Uri{s)Ci, (2.4) 

1=1 

where 

Uri{s) = Ui{s) - Ui,{s)l'i{s). (2.5) 

We consider the following boundary conditions on ^i: 

u^(s) = 0, se Si, (2.6) 
F^s) = -x{F,{s), \u^s)fX{s), s e Si. (2.7) 

Here % is the function of slip that depends on the normal component of the surface force F^, 

and on the square of the module of the tangential velocity u'^. 

Formula (2.7) is a generalization of Navier's condition of slip in which X is a positive 
constant, the nonlinear modiflcation of Navier's condition of slip in which x is a function 
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of and Coulomb's law of friction in which x = oo at |F^| < cilF'^j and x = c at 
\F'^\ = cilF'^l, c,ci are positive constants. 

Wc note that problems on flow of nonlinear viscous fluids in which % is a function of l-u"^! 
were investigated in [8]. 

The function x accepts positive values, x does not depend of Fi, at Fi, > 0, and it rises 
as decreases. The sign minus in (2.7) designates that the velocity of slip of the fluid is in 
opposition to the tangential surface force, i.e. the frictional force is in opposition to the direc- 
tion of motion, and the module of the slip velocity is equal to |F'^(s)|(x(F,y(s), \u'^ {s)\'^))~^ . 

In the special case that x(yi,y2) = oo for an arbitrary (2/1,2/2) £ ^ x formulas 



= 0, i.e. the fluid adheres to the hard boundary, and in the case of 



(2.6), (2.7) imply u 

X{yi^y2) = for an arbitrary (2/1,2/2) G M x M+, (2.7) yields F 



= 0, i.e. the frictional 

Si 

force is equal to zero. The relation (2.6) designates that the fluid does not flow through the 
hard wall Si. 

For the constitutive equation (1.1) the components Fi of the surface force F = (Fi, . . . , F„) 
are defined by 

Fi = [-pSij + 2(p{I{u), \E\,fj,{u,E))eij{u)]uj i = l,...,n, (2.8) 

Si 

and the normal component F^, of the surface force is determined by (2.2). In (2.8) and below 
the Einstein convention on summation over repeated index is applied. 
Let P be an operator of regularization given by 

Pv{x)= / io{\x - x'\)v{x') dx' , xeU, (2.9) 

where 

a; G C°° (M+), suppo; G [0,a], u{z) > 0, ze R+, 

u}{\x\) dx = 1, a is a small positive constant. (2-10) 



/ 



Here, we assume that the function v is extended to W\ 

We denote by Fru(j>,u) the normal component of the surface force calculated by the 
regularized functions of pressure p and velocity u. According to (2.2) and (2.8), the function 
Fri,{p,u) is defined as follows: 

Fr^{p,u) = [-Pp + 2^{I{Pu),\E\,^i{Pu,E))eij{Pu)uiUj] . (2.11) 

Si 

We change the function Fi, in (2.7) for the function Fj.i/{p, u). Then, we obtain the following 
boundary condition: 

F^s) = -x{Fru{p,u){s),\u^is)\^Xis), seSi. (2.12) 

From the physical point of view, (2.12) denotes that the model is not local, the velocity 
of slip at a point s & Si depends on the averaged normal surface force Fru{p,u) which 
in its turn is defined by the values of pressure and the derivatives of the velocity at points 
belonging to some small vicinity of the point s. This is natural from the physical view-point. 

Such nonlocal approach is also connected with the fact that the velocity of slip depends 
on the surface roughness which is not a local characteristic. 
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Taking (2.2), (2.3) and (2.4) into account, we represent (2.12) in the following form: 

n n 
Fi- {^FkVkjl'i = -x{Frv{p,u),^ul^Uri OTL Si, 1 = 1, . . . ,71. 

k=i k=i (2-13) 

Finally, we obtain by (2.8) and (2.13) the following boundary condition of slip: 

n 

2ip{I{u),\E\, n{u, E))[eij {u)iyj - ekj{u)vjVkVi] = -x{Fr„{p,u), '^ul^)uri on Si, 

k=l 

i = l,...,n, (2.14) 

where Fri,{p,u) is defined by (2.11). 

3. Governing equations and assumptions. 

We consider stationary flow problem under the Stokes approximation, i.e. we ignore 
inertial forces which are assumed to be small as compared with the internal forces caused 
by the viscous stresses. Then, the motion equations take the following form: 

dv d 

■^~2—[ip{I{u),\E\,ii{u,E))eij{u)] = Kiiun, i = l,...,n, (3.1) 

where Ki are the components of the volume force vector K. 
The velocity function u meets the incompressibility condition 

div« = y^ = in a (3.2) 

OXi 

i=l 

We assume that Q is a bounded domain in M", ra = 2 or 3 with a Lipschitz continuous 
boundary S. Suppose that Si and 52 are open non-empty subsets of S such that Sir]S2 = 0, 
and 5i U ^2 = S. 

We consider mixed boundary conditions for the functions u, p. Wherein, the terms of slip 
(2.6), (2.14) are specified on Si and surface forces are given on 5*2, i.e. 

[-pSij + 2(p{I{u),\E\,iJ,{u,E))eij{u)]iyj = Fi, i = l,...,n. (3.3) 

We assume that 

(CI): if : {yi,y2:yz) 95(2/1) 2/2, ys) is a function continuous in x [0, 1], and for an 
arbitrarily fixed (2/2,2/3) G M+ x [0,1] the function ip{.,y2,y3) ■ yi ¥'(2/1,2/2,2/3) is 
continuously differentiable in R+, and the following inequalities hold: 

02 > (^(2/1,2/2,2/3) > oi (3.4) 
¥^(2/1,2/2,2/3) + 2 — (2/1,2/2,2/3)2/1 > «3 (3.5) 

oyi 

^ (2/1,2/2,2/3) yi < 04, (3.6) 
oyi 

where a^, 1 < i < 4, are positive numbers. 

Inequality (3.4) indicates that the viscosity is bounded from below and from above by 
positive constants. The inequality (3.5) implies that for fixed values of \E\ and fi{u,E) the 
derivative of the function I{v) — > G{v) is positive, where G{v) is the second invariant of the 
stress deviator 

G{v) = A[<p{I{v),\E\, fi{u,E))]''l{v). 
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This means that in the case of simple shear flow the shear stress increases with increasing 

dip 
dyi 



shear rate. (3.6) is a restriction on ^ for large values of yi. These inequalities are natural 



from the physical point of view. 

Relative to the function of slip x, we assume that the following conditions are satisfied: 
(C2): X '■ (yi) 1/2) x(yi) 2/2) is a function continuous in M x M_|_, and for an arbitrarily 
fixed yi G M, the function xiVi^ ■) '■ y2 ^ x{yi-iy2) is continuously differentiable in 
R+, and the following inequalities hold: 

^'2 > x(yi,2/2) > (3.7) 

X(yi,y2) + 2— (yi,y2)?/2 > ^3, (3.8) 
Qx 

^(2/1,2/2) 2/2 < &4, (3.9) 

where (2/1,2/2) G R x M_|_, and bi, 1 < i < 4, are positive numbers. Inequalities (3.7)- 
(3.9) are analogous to the ones of (3.4)-(3.6). Inequality (3.7) means that the function 
of slip is bounded from below and from above by positive constants. (3.8) implies that 

for fixed value of yi, i.e. the value of Fruip,u){s), s £ Si, the derivative of the function 

1 

\u'^\ \F'^\ = x{Fruip,u), |tt'^p)|n'^| is not less than 63. Indeed, denoting z = y2 , we obtain 

Q-^ =x{yi,z)+2—{yi,z)z >h, 

that is the frictional force increases as the velocity of slip increases. 

The inequality (3.9) is a restriction on for large values of 2/2- The inequalities (3.7)— 
(3.9) arc natural from the physical viewpoint. 

We suppose also that 

K = {Ki,...,Kn)eL2{nT, F = (Fi,...,F„)GL2(52)". (3.10) 
4. Boundary value problem. 

We study a problem of searching for a pair of functions {u,p) which satisfy the motion 
equation (3.1), the condition of incompressibility (3.2) and the boundary conditions (2.6), 
(2.14) and (3.3). 

Consider the following spaces 

Z = {v\veH\nr, vAs,=^}, (4.1) 
W = {v\veZ, divv = 0}. (4.2) 

Lemma 4.1. Let he a hounded domain in R", n = 2 or 3 with a Lipschitz continuous 
boundary S, and let Si he an open nonempty subset of S. Then the expression 

||^;||z= ( / I{v)dx+ / Y.^'rhdsY (4.3) 

^ J Si ^^-^ ' 

defines a norm in Z and W being equivalent to the norm of H^{^)^. 
For a proof see in [7], Section 1.7. 

Everywhere below we use the following notations: If F is a normed space, we denote by 
Y* the dual of Y, and by (/, h) the duality between Y* and Y, where f & Y* , h ^ Y . In 
particular, if / G L2{^) or / G L2(f^)", then (/, /i) is the scalar product in L2(fi) or in 
L2(0)", respectively. The sign ^ denotes weak convergence in a Banach space. 
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Denote by B the operator of divergence, i.e. 

Bu = divu. (4.4) 

It is obvious that B is a hnear continuous mapping of Z into L2{0,), i.e. B G C{Z, L2{'^)). 
We denote by B* the adjoint to B operator. 

We introduce operators M : Z ^ Z* and A : Z x L2{^) ^ Z* as follows: 

{M{u),h) = 2 (p{I{u),\E\,n{u,E))eij{u)eij{h)dx, u,heZ, (4.5) 
Jn 

„ n 

{A{u,p),h)= x(Fru{p,u),y2^lk)'^rihTids, {u,p) e Z X L2{n), heZ. 

•Jsi ^ k=i ' (4.6) 

Consider the problem: find a pair (u^p) such that 

u^Z, peL2in), (4.7) 

iM{u),h) + {A{u,p),h)-{B*p,h) = {K + F,h), heZ, (4.8) 

(Bu,5) = 0, qeL2{n). (4.9) 

Here we use the notations 

iK,h)= I K.h^dx, {F,h)= [ Fihids. (4.10) 

Jo. JS2 

A solution of the problem (4.7)-(4.9) will be called a generalized solution of the problem 
(3.1), (3.2), (3.3), (2.6) and (2.14). Indeed, by use of Green's formula it can be seen that, if 
(n,p) is a solution of the problem (3.1), (3.2), (3.3), (2.6) and (2.14), then (n,p) is a solution 
of the problem (4.7)-(4.9). On the contrary, if (n,p) is a solution of the problem (4.7)— (4.9), 
then (n,p) is a solution of the problem (3.1), (3.2), (3.3), (2.6) and (2.14) in the sense of 
distributions. 

Theorem 4.1. Let Q, he a hounded domain in M", n = 2 or 3, with a Lipschitz continuous 
houndary S, and suppose that the conditions (CI), (C2) and (3.10) are satisfied. Then there 
exists a solution of the problem (4.7) -(4.9). 

5. Auxiliary results. 
We consider four functions vi, V2, vs, such that 

vi G L2{0,), vi{x) > a.e. in fi, 
V2 G L^{n), V2{x) e [0, 1] a.e. in fl, vs G L2{n), G H^ip.^. (5.1) 
We set V = {vi,V2,Vi, V4) and define the operator My : Z ^ Z* as follows: 

(My(u),e) = 2 / 'f{I{u),vi,V2)eij{u)eij{e)dx + j x {Fru{v3, "^^4), ^rfc) '^ri ^ri ds^ 
JQ J Si k=l 

u,eeZ. (5.2) 

Lemma 5.1. Suppose that the conditions (CI), (C2) and (5.1) are satisfied. Then 

{M^{u) — My{w),u — w) > niWu — w\\'^z^ u,w E Z, (5-3) 
\\My{u) - My{w)\\z* < H2\\u- w\\z, u,w e Z, (5.4) 



where 

fii = min(2ai, 203,61,63), 

fi2 = 202 + 4a4 + 62 + 264. (5.5) 
Proof. We present the operator in the form 

= Ml + M2, (5.6) 

{Mi{u),e)=2 (p{I{u),vi,V2)eij{u)eij{e)dx, (5.7) 
Jn 

r. n 

Urierids, u,eE. Z. (5-8) 

•^■^1 k=i 
Let u, w be arbitrary functions in Z and 

h = u — w. (5-9) 

We introduce the function 7 as follows: 

7(t) = / ip{I{w + th), vi,V2)eij{w + th)eij{e) dx, t e[0,l], e e Z. 
Jn 



(5.10) 



It is obvious that 



7(1) - 7(0) = ^(Mi(u) - Miiw),e). (5.11) 

By using the theorem on the differentiability of a function represented as an integral, we 
conclude that 7 is differentiable at any point t G (0, 1). Therefore 

7(1) = 7(0) + ^(0, ^G(0,1), (5.12) 



where 



dj f 

-^iO = JJ_ip{I{w + ^h),vi,V2)eij{h)eij{e) 



dip 

+2— (/(u; + ^/i),ui,?;2)£fem(^i^ + ?^)£fcm(^)£ij(^i^ + ?^)£y(e)]c?a;. (5.13) 

Taking note of the inequality 

\ekm{w + £,h)ekm{h)\<I{w + ih)'2l{h)'2, (5.14) 
and (3.4), (3.6), (5.11)-(5.13), we obtain 

\\Mi{u) - Mi{w)\\z* < (2a2 + 4a4)( J I{u-w)dxy 

< (2a2 +4a4)||n-u;||z. (5.15) 

Define the function g as follows: 

( X /^(«'^i(^)'^2(a;)), ^ia,vi{x),V2ix)) <0, 
^^"'^^ \0, if ^(a,^;i(x),^;2(x))>0, 

where a G M+, x E 0,. 



Then, taking e = h in (5.13) and applying (3.4), (3.5), (5.9) and (5.14), we obtain 

dt 



"(6= / Hl{w + ^h),Vl,V2)I{h) 

Jn 



dip f 
+2- — (/(ui + ^/i),t)i,f2)(eij(iL' + ^/i)eij(/i))^]da:; > min(ai,a3) / I{u-w)dx 

oyi Jn (5.16) 
(5.11), (5.12) and (5.16) imply 

{Mi{u) - Mi{w),u-w) >2mm{ai,a3) / I{u-w)dx. (5.17) 
We introduce the function 71 as follows: 

» n 

7l(*) = / x(^Fruiv3,V4),'^{Wrk + thrkf^iWri + thri)eri ds, t G [0,1], e E Z, 

•^•^1 k=l 

where h is defined by (5.9). 

By analogy with the foregoing, we obtain 



(M2(«) - M2{w),e) = 71 (1) - 71 (0) = ^(6) 



dr/i 
dt 



X (u^rfc + iihTk)hTk{wTi + iihri)e-Ti ds, G (0, 1), (5.18) 
and (3.7)-(3.9) imply 

an 1 
_ ^{Uri - Wrif dsy 
-'i 1=1 

< (62 + 264)||n-w;||z, (5.19) 

„ n 

(M2(tx) - M2(i(;), u - u;) > min(6i, 63) / ^(u^fe - "i^rfc)^ t^s- (5-20) 

"''S'l fc=i 

Taking (5.15), (5.17), (5.19) and (5.20) into account, we obtain (5.3)-(5.5). ■ 
Let a be a positive number. Define the operator : Z ^ Z* as follows: 

f 1 " 
{Aaiu),h)= x(^rv( Bu,u),^ulf?jurihTids u,heZ. (5.21) 

•^■^1 " k=l 

Consider the problem: find a function Ua satisfying 

ua G Z, (5.22) 

{M{ua),h) + {Aa{ua),h) + -{Bua,Bh) = {K + F,h), heZ. (5.23) 

a 

The problem (5.22), (5.23) is an approximation of the problem (4.7)-(4.9) in which the 
function of pressure p is replaced by the function —a~^ divu; in this case we do not assume 
that div-Ua = 0. 

Theorem 5.1. Let Q be a bounded in M", n = 2 or 3, with a Lipschitz continuous boundary 
S. Suppose that the conditions (CI), (C2) and (3.10) are satisfied. Then for an arbitrary 
a > 0, there exists a solution of the problem (5.22), (5.23). 
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Proof. Let {Zm}m=i ^'^ ^ sequence of finite dimensional subspaces in Z such that 

hm inf \\v - h\\z = 0, v e Z, (5.24) 
Zm C Zm+1, men. (5.25) 
We seek an approximate solution of the problem (5.22), (5.23) in the form 

^am € ^rri! (-^(^am))^) + {-^ai^am 

),h) + -{Buam,Bh) = (K + F,h), h e Zm- 

" (5.26) 

By (3.4), (3.7), (3.10) and (4.3), we obtain 

y{h) = iM{h), h) + (A„(/i), h) + ^{Bh, Bh) -{K + F, h) 

>2ai f I{h)dx + bi [ y2'^rids-\\K + F\\z*\\h\\z>tii\\hfz-\\K + F\\z*\\h\\z, 

-^Si i=i (5.27) 

where /xi = min(2ai, bi). 

Therefore, y{h) > for \\h\\z > r = \\K + FWz* /i^^- 

From the corollary of Brouwer's fixed point theorem (cf.[2]), it follows that there exists a 
solution of (5.26) with 

\\Uam\\z<r, \\M{Uam) + Aa{Uam)\\z* < C, m G N, 

where the second inequality follows from (3.4) and (3.7). Therefore, we can extract a 
subsequence {uq,^}^^^ such that 

Z, (5.28) 

Uar) Ua in L2(f^)"' and a.e. in J7, (5.29) 

M(m«^) + A^iuarj) ^ 6* in Z*. (5.30) 

Let 770 be a fixed positive integer and h G Z^^. Observing (5.29), (5.30), we pass to the limit 
in (5.26) with m replaced by rj, and obtain 

(e + -B* Bua,h) = (K + F,h), heZrjo. (5.31) 
a 

Since 770 is an arbitrary positive integer, by (5.24), we obtain 

e + -B*Bua = K + F inZ*. (5.32) 
a 

We present the operators M and in the form 

M{u) = M{u,u), Aa{u) = Ao,{u,u), (5.33) 

where the operators {u,v) — M{u,v) and {u,v) — Aa{u,v) are mappings of Z x Z into Z* 
according to 

{M{u,v),h) =2 [ <p{I{v),\E\,fi{u,E))eij{v)eij{h)dx, (5.34) 
Jn 

C 1 " 

{Aa{u,v),h) = J x{Frv{--Bu,u),^^vl^Vrihrids. (5.35) 



k=l 
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Denote 



-M{Uari,v) ~ Aa{Uarj,v),UaT, - v), V E Z. (5.36) 

By Lemma 5.1, see (5.3), we obtain 

Xniv) > 0, r? G N, v e Z. (5.37) 

We have 

\\M{Uar„v)- M{Ua,v)\\z* 

< 2{ jWiv), \Eln{uar,,E)) - ^{I{v), \E\,n{ua,E))]^ I{v) dxy. (5.38) 

(CI), (5.29), (5.38), and the Lebesgue theorem give 

M{uar,,v) ^ M{ua,v) in Z* . (5.39) 

Likewise, we obtain 

AaiUar],v) AaiUa,v) in Z* . (5.40) 

(5.26), (5.28) and (5.33) yield 

" (5.41) 

By (5.30) and (5.32), we obtain 

linLL[{M{Uar,,Uar,) + Aa{Ua,j,Ua,j),v)] + -{BUa,Bv) = {K + F,v), V E Z. 

" (5.42) 
Observing (5.39)-(5.42), we pass to the limit in (5.36). Then by (5.37), we find 

(K + F-M(ua,v)-AJua,v)--B*Bv,Ua-v)>0, veZ. (5.43) 

a 

We choose v = — ^h, j > 0, h e Z, and consider 7 — 0. Then, Lemma 5.1, see (5.2), 
(5.4), (5.34) and (5.35), implies 

{K + F- M{ua) - Aa{ua) -^B* Bua, h) > 0. (5.44) 

This inequality holds for any h e Z. Replacing hhy —h shows that equality holds true in 
(5.44). Therefore, n„ is a solution of the problem (5.22), (5.23). ■ 

We will use also the following lemma: 

Lemma 5.2. Let be a bounded domain in M", n = 2 or 3, with a Lipschitz continuous 
boundary S, and let the operator B G jO{Z,L2{^)) be defined by (4.4). Then, the inf-sup 
condition 

inf sup „ i/^i'^'f^ > ^1 > (5.45) 

MGi2(n) v(zz \\v\\z ||M||L2(n) 

holds true. The operator B is an isomorphism from W-^ onto L2{^), where W-^ is orthogonal 
complement of W in Z, and the operator B* that is adjoint to B, is an isomorphism from 
L2{p) onto the polar set 

W^o = {/ G Z*, (/, u) = 0, u G W}. (5.46) 
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Moreover, 



B ^\\c{L2{n),w-^) ^ ^> (5-47) 
1 



\\iB*r'\\c(WOMm<^- (5.48) 



Z*\\Uai\\z- 

(6.4) 



For a proof see in [7], Section 6.1.2. 

6. Proof of Theorem 4.1. 

Let {ui} be a sequence of positive numbers such that hmaj = 0. Consider the problem: 
given ccj, find -Uq. satisfying 

e Z, (6.1) 

(M(«„,),/i) + (A„,(t/aJ,/i) + ar^(^«",,S/i) = (i^ + i^,/i), heZ. (6.2) 

The existence of a solution of the problem (6.1), (6.2) follows from the Theorem 5.1. 
Taking h = in (6.2), we obtain 

{M{Uai),Uai) + {Aai{Ua,),Ua,) + a~^\\BUa,\\l^(^^) < + i^||z* ll^ai lU- (6.3) 

(3.4), (3.7) and (6.3) imply 

2ai f I{ue,,)dx + bi I ^ul.^J^ds + ar^\\Bua^\\l^^a)<\\K + F 
Jo, J Si 

It follows from here and (4.3) that 

II^qJIz < ci, (6.5) 
_ 1 

«i 'l|-BWaJ|L2(n) < C2. (6.6) 

Therefore, a subsequence {ua^ } can be extracted from the sequence {w^ . } such that 

Uam ^ in 2", (6.7) 

Uam ^ in L2{^i) and a.e. in J7, (6-8) 

Uamlsi ^ u\si in L2(5'i) and a.e. in 5i, (6.9) 

Bua^^O mL2{n). (6.10) 

(6.2) yields 

a-}B*Bu,,„^ = K + F-M{u^J-Aa^{uaJ inZ*. (6.11) 

By virtue of (3.4), (3.7), (6.5) and (6.9) the right-hand side of (6.11) is bounded in Z*. 
Therefore, Lemma 5.2, yields 

(^m\\BUaJ\L2{n) < C3, (6.12) 

and we can consider that 

a:;^Bua^^p inL2(ri). (6.13) 

By analogy with the proof of Theorem 5.1, wc pass to the limit in (6.2) using (6.7)-(6.10) 
and (6.13). As a result, we obtain that the pair {u,p) is a solution of the problem (4.7)-(4.9). 
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Remark. Suppose that the condition of shp has the following form: 

2ip{I{u),\E\,ii{u,E)) eij(u)vj - ekj{u)vjVk^i 

n 

= -X{FrAPi ^{PUrkf)Uri On ^i, i = l,...,n, (6.14) 
k=l 

compare with (2.14). Here P is the operator of regularization defined by (2.9), (2.10). 

In this case, the function x in the operator A in (4.6) is defined by just the same expression 
as in (6.14). 

The condition of slip (6.14) is reasonable from the physical point of view and under such 
condition the Theorem 5.1 remains true without the restrictions (3.8) and (3.9). 

7. Galerkin method for the problem (4.7)-(4.9). 
Let {Njn} be a sequence of finite-dimensional subspaces in L2{^) such that: 

lim inf \\w - y\\L^/Q) = 0, weL2{^). (7.1) 

NraCNrn+1, m G N, (7.2) 

inf sup „ l^irY^ > /? > 0, m G N, (7.3) 
where the operators B„i € C{Zm,N^) are defined as follows: 

{BmV,iJ,)= / iJ,divvdx V e Zm, n^Nm- (7.4) 

Let B*^ G £{Nm,Z^) be the adjoint operator of Bj^ with {B^v^ji) = (v,5^/x) for all 

V G Z„i and all ^ G N^. 

We introduce the spaces Wm and by: 

W^ = {v£ Zra, [BmV, /i) = 0, ;U G iY^}, (7.5) 

Wl = {q^Zl,, (g,i;) = 0, v<-WrA. (7.6) 

Lemma 7.1. Let {Zm}m=i> {^m}m=i sequences of finite- dimensional subspaces in Z 
and L2{fl) such that (7.3) holds true. Then the operator B^ is an isomorphism from 
onto W^, and the operator Bjn is an isomorphism from onto N^, where is an 
orthogonal complement ofWm in Z^. Moreover, 

ll(^m)"'ll/:(w^o,;v^)<;^, \\^m\\t{N*^,w^)<-^. ^eN. (7.7) 
For a proof see in [7], Section 6.1.2. 

We seek an approximate solution of the problem (4.7)~(4.9) of the form: 

iUm,Pm) G X Nm, (7.8) 

{M{um),h) + iA{um,Pm),h)-iB*^Pm,h) = {K + F,h) heZm, (7.9) 
{BmUm,q)=0, qeNm, (7.10) 

Theorem 7.1. Let U be a bounded domain in M", n = 2 or 3, with a Lipschitz continuous 
boundary S, and suppose that the conditions (CI), (C2) and (3.10) are satisfied. Let also 
{Zrn}m=i ^'^'^ {^rn}m=i Sequences of finite- dimensional subspaces in Z and L2{^) which 
satisfy the conditions (5.24), (5.25), (7.1)-(7.3). Then for an arbitrary m G N, there exists 



14 

a solution of the problem (7.8)-(7.10), and a subsequence {tffejPfc} can be extracted from the 
sequence {um,Pm\ such that: 

Uk ^ u in Z, (7-11) 
Pk^p inL2{n), (7.12) 

where {u,p) is a solution of the problem (4.7)-(4.9). 

Proof. 1) We determine a mapping Ai : L2{^1) x Z x Z — as follows: 

„ n 
{Ai{lJ,,V,w),h) = / xi^ri^ifJ-yV),^ wlk)Wrihrids, 
"'■S'l k=l 

fi<^L2in), {v,w,h)eZ^. (7.13) 

Consider the following the problem: given a pair {vm, IJ-m) G -^m x ^m, find {vm,fi'm) 
satisfying 

X Nm, (7.14) 
(M(t;,„,?)^),/i) + (Ai (//,„, i;^, {),„), /i) - {B^fim,h) = {K + F,h), he Zm, 

(7.15) 

{Bm,Vm,q)=0, qeNm. (7.16) 
It follows from (7.5) and (7.14)-(7.16) that Vm is a solution of the problem: 

Vm G Wm, (7.17) 

(M(i;^,«^),/i) + (Ai(/x^,?;^,«^),/i) = (ii: + F,/i) /i G W"„. (7.18) 
By Lemma 5.1, the operator 

M{Vm, ■) +Ai{lJ,m,Vm, ■) : V ^ M{Vm,v) + Ai{lJ.m,Vm,v) 

is a mapping of Z into Z* that is strictly monotone, coercive and continuous. Therefore, 
there exists a unique solution of the problem (7.17), (7.18), and by Lemma 7.1, there exists 
a unique function flm G Nm such that: 

B^P-m = M{Vm,Vm) + Ai{lJ,m,Vm,Vm) - K - F in Z^. (7.19) 

In this case the pair (vm, P'm) is a unique solution of the problem (7.14)— (7.16). 
We take /i = ■0^ in (7.18). Then by (3.4) and (3.7), we obtain 

\\Vm\\z<Ci, {Vm, IJ-m) & Zm X Nm, m G N. (7.20) 

Lemma 7.1 and (7.19) imply 

||/im||L2(n) < C2, {Vm,IJ'm) & Zm X Nm, m G N. (7.21) 

For m G N, we introduce a mapping '. X Nm — ^ ^ follows: ('^^m? j-^rri) ^ 

Zm X -/Vm, Bmivm,IJ'm) = (*m,/im), where (i)m,Am) is the Solution of the problem (7.14)- 
(7.16). 

Let {gk,ak} e Z^ x and gu ^ g, ak ^ a. By using (CI), (C2), (7.20), and (7.21), 
one can verify that Bm{gki ctk) ^ I3m{g, ct). Hence, Bm is a continuous mapping of Zm x Nm 
into itself. 

Moreover, (7.20) and (7.21) yield that the mapping Bm maps a compact convex set 

dm = {{v,IJ-) e ZmX Nm, \\v\\z < Cl, ||/x||l2(0) < C2} 



15 



into itself. Therefore, the Schauder principle implies that there exists a pair {um,Pm) £ 



X Nm such that: 

Sm{Um,Pm) = {Um,Pm)- (7.22) 

In addition, the pair {um,Pm) is a solution of the problem (7.8)-(7.10) for any m, and we 
have 

\\um\\z < ci, ||Pm||L2(n) < ^2, m € N. (7.23) 

Hence, a subsequence {uk,Pk} can be extracted from the sequence {um,Pm} such that: 

Uk no in Z, (7.24) 

Uk — ^ uo in L2{^) and a.e. in J7, (7-25) 

Pk^Po inL2{0,), (7.26) 

Fru{Pk,Uk) ^ FrAPo,uo) in Loo(S'i), (7.27) 

M{uk) + A{uk,Pk)^@ inZ*. (7.28) 

Let ko be a fixed positive integer and let h G Z^,,, q G A^'^,,. By (7.24), (7.26), (7.28), we 
pass to the limit in (7.9), (7.10) with m changed by k, which gives 

{e-B*po,h) = {K + F,h), heZk,, (7.29) 

/ qdivuodx = 0, qeNk^. (7.30) 
Jn 

Since ko is an arbitrary positive integer, we obtain by (5.24), (7.1), (7.29) and (7.30) that 

e-B*po = K + F, inZ* (7.31) 

divno = 0. (7.32) 
We determine a mapping Jj. : Z — > Z* as follows: 

{Jk{v),h) = {M{uk,v)+Ai{pk,Uk,v),h), k = 0,1,2,.... (7.33) 
It follows from (5.34) and (7.13) that 

Jkiuk) = M{uk) + A{uk,Pk), (7.34) 

and Lemma 5.1 gives 

{Jk{uk)- Jk{v),Uk-v)>0, V e Z, fc = 0,1,2, (7.35) 

(7.25), (7.27) and the Lebesgue theorem imply 

lim{Jk{v),Uk) = {Jo{v),uo), 

lim( Jjfc {v) ,v) = Jo {v) , v) , (7.36) 
By (7.28), (7.31) and (7.34), we have 

\im{Jk{uk),v) - {B*po,v) = {K + F,v). (7.37) 

Taking into account that 

{BkPk,Uk) = {Pk,Buk) = 0, 
we get from (7.9), (7.24) and (7.34) that 

(Jkiuk), Uk) = {K + F, Uk) ^{K + F, uo). (7.38) 
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Upon (7.36)-(7.38), we pass to the limit in (7.35), which gives 

{K + F-Jo{v) + B*po,uo-v)>0, veZ. (7.39) 

Take here v = uq — ^h, ^ > 0, h e Z, and let ^ tend to zero. By Lemma 5.1, see (5.4), we 
get 

{K + F- Joiuo)+B*po,h)>0, heZ. (7.40) 

Therefore, the pair u = uq, p = po is a solution of the problem (4.7)-(4.9). 

2) We will show that the solution of the problem (7.8)-(7.10) converge to the solution of 
(4.7)-(4.9) strongly. 

Let 

Yk = {Jk{uk) - M'^o), Uk - uo). (7.41) 

Obviously 

Yk = (Jkiuk) - Jk{uo),Uk - Uq) + (Jfc(no) - Joiuo),Uk - Uq). (7.42) 

Upon (7.24), (7.36)-(7.38) and (7.41) limYk = 0. By (7.24), (7.36) the second addend in 
(7.42) tends also to zero, and so 

\im{Jk{uk) - Jk{uo),Uk - Uq) = 0. (7.43) 

Observing (7.33), (7.43) and Lemma 5.1 sec (5.3), wc obtain (7.11). 

We take /i G in (4.8) and (7.9) and subtract (4.8) from (7.9). Then, we obtain 

{B*{pk - /x), h) = {M{uk) + A{uk,Pk) - M{uq) - A{uQ,po), h) + {B*{pq - n), h), 

heZk, /X G Nk. (7.44) 

This equality together with (7.3) yields 

(B^ipk m) ^) 

\\Pk - t^Wh^iQ) < sup — ■ — < Ak + c\\po - /i||L2(n), M e Nk, 

heZk PlMz (^745) 

where 

Ak= \\M{uk) + A{uk,Pk) - M{uo) - A{uo,po)\\z*. (7.46) 

Hence 

WPo-Pkh^n < inf (Ibo - MllL2(n)) + \\Pk - /^Ikam)) 

<r'A + (c+l) inf ||po-/x|U,(n). (7.47) 

By (7.11) and (7.26), we obtain lim^^ = 0, and (7.1) implies 



lim inf ||po - /^IlLzfQ) = 0> 



so that (7.12) follows from (7.47). 
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